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Abstract 

We show that recent observations of fractal dimensions in the //-space of iV-body 
Hamiltonian systems with long-range interactions are due to finite N and finite 
resolution effects. We provide strong numerical evidence that, in the continuum 
(Vlasov) limit, a set which initially is not a fractal (e.g. a line in 2D) remains 
such for all finite times. We perform this analysis for the Hamiltonian Mean Field 
(HMF) model, which describes the motion of a system of N fully coupled rotors. 
The analysis can be indirectly confirmed by studying the evolution of a large set of 
initial points for the Chirikov standard map. 
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It has been recently claimed [1,2] that structures characterized by a non integer 
fractal dimension [3] form spontaneously in the /i-space of some iV-body Hamiltonian 
systems, even when starting from a set of initial conditions (points in /i-space) that 
does not show any fractality. 
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The model analysed in Ref. [1] is the one-dimensional self-gravitating system [4] and 
the initial condition is a line in /i-space, corresponding to a vanishing virial ratio 
(zero velocity for all mass sheets). The model studied in Ref. [2] describes the motion 
of a fully coupled system of N rotors and goes under the name of Hamiltonian Mean 
Field (HMF) model [5] ; the initial condition is such that all rotors are fixed at some 
given angle (conventionally zero) and angular velocities are uniformly spread in an 
interval, symmetric around zero. 

There is a strong analogy between the two results since: i) both models possess 
a continuum mean-field limit (N — > oo at fixed volume) for which the single par- 
ticle distribution function (SPDF) f(q,p,t) evolves according to a Vlasov-Poisson 
system [6], ii) both initial states are lines in /i-space when the continuum limit is 
taken. 

In this short note we analyse the HMF model, considering the same initial state 
as the one studied in Ref. [2]. We give compelling numerical evidence that, as the 
number of rotors grows and as the resolution increases, the collection of initial points 
lies on a set of dimension one in the /i-space (which is two-dimensional for this model) 
at all finite times. 

We will explain that this is a somewhat trivial consequence of the fact that the 
SPDF obeys a Vlasov equation in the mean-field limit. Moreover, we will show that 
apparent fractal dimensions can be indeed observed at low resolutions, justifying the 
observations of Ref. [2]. For comparison, we will show what happens for the standard 
Chirikov's map [7], which shares the symplectic property with the HMF model, for 
which we can push the analysis to a larger number of initial points. 

The Hamiltonian of the HMF model is 

I jv 1 N 

Hhmf = - $>? + ^ E t 1 - cos fe - Qk)}, (1) 

j=l j,k=l 

where G [— 7r,7r[ is the position (angle) of the 2-th particle on a circle and Pi the 
corresponding conjugate variable. This system can be seen as representing particles 
moving on a unit circle interacting via an infinite range attractive cosine potential, 
or as classical XY rotors with infinite range ferromagnetic couplings (for more details 
see Ref. [8]). The magnetization, defined as 

-> 1 N 

M{t) = (M x , M y ) = - Y / (cosq J (t),smq j (t)), (2) 

iV 3=1 

is the main observable that characterizes the dynamical and thermodynamic state 
of the system. 

In the continuum limit, that is keeping the volume (here the interval [— 7T, 7r[) and the 
energy per particle fixed as the number of particles N — > oo, the dynamics governed 
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by Hamiltonian (1) is described by a Vlasov equation. Indeed, the state of the finite 
N system can be described by a single particle time-dependent distribution function 



1 N 

fd (q,P, t) = — s (<? - Qj (*) ,P - Pj (*)) > ( 3 ) 

where 5 is the Dirac function. When N is large, it is natural to approximate the 
discrete density fd by a continuous SPDF / (q,p,t). Using this distribution, one can 
rewrite the two components of the magnetization M as 



M x [/] = J f(q,p,t) cosg dqdp , (4) 
M y [f]=J f(q,p,t) sinqdqdp . (5) 

Within this approximation the potential that affects all the particles is 

V{q)\f\ = l-M x [f\ cos q-M y [f] sing . (6) 
This potential enters the expression of the Vlasov equation 

df df dV df 

m +p d- q -d q - [f] dp- = ° ' (7) 

which governs the spatio-temporal evolution of the SPDF /. 

The initial condition considered in [2] is qt = 0, Vz and pi random i.i.d. uniformly 
distributed in the interval [—p, p] ■ As usual, we call this class of initial conditions 
water bags (WB). The value of p is chosen such that the energy per particle is 
U = Hhmf/N = 0.69 and N = 10000, but these data will be varied in this paper. 
For such an energy this initial state is Vlasov unstable [9], and hence the SPDF 
evolves in time (see fig. 3 in Ref. [2]). 

However, since initially the state is a line in the two-dimensional /x-space (q,p), it 
will remain a line for all finite times, although intricately stretched and folded. This 
is a trivial property related to the fact that Vlasov equation defines a flow $ t for 
the particles (points x = (q, p) in the /x-space continuum) 

x(t) = * t (x(0)) , (8) 

which is a diffeomorfism. Moreover, since the Jacobian J = <9x(£)/<9x(0) is a sym- 
plectic matrix, areas modify their shape but are conserved. 

The reason why the set does not appear to be a line is due to the finite particle 
number in the simulation. The line is stretched and folded by the two main phase- 
space mechanisms: the formation of whorls around elliptic fixed points and of tendrils 
in the homoclinic tangle [10]. Moreover, local hyperbolicity in phase-space produces 
a dispersion of the initial points. As a consequence, the set of initial points, appears 
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by visual inspection, in the course of time, to be partly distributed along lines and 
partly dispersed onto the two-dimensional /i-space (see again fig. 3 in Ref. [2] or fig. 
2 in Ref. [1]) This is why the authors of Refs. [1,2] have used fractal dimension to 
chacterize this set, more specifically box dimension [1] and correlation dimension [2]. 

The number of boxes iV&(e) of linear size e which cover the set of N = 10 6 points at 
time t — 75 when starting the HMF model in a WB initial condition at U = 0.69 
is shown as a function of resolution 1/e in log- log scale in fig. 1. These data are 
not qualitatively different from the ones in fig. 3 of Ref. [1] and one might then fit 
the central part of them with the Ansatz Nf, = e~ D ° (since the maximal box size 
which contains all the set is normalized to one), obtaining D Q ~ 1.930 ± 0.002 and 
concluding that the set is a fractal, with dimension D close to two, but definitely 
distinct from two within numerical accuracy. This is a consequence of the set being a 
mixture of a set with one- dimensional features (the initial line) and two-dimensional 
ones caused by chaotic spread. 



Fig. 1. Number of occupied boxes of linear size e versus resolution 1/e in log-log scale for the 
HMF model in a WB initial condition with U = 0.69 and N = 10 6 (<>)• The short-dashed 
line superposed to the data is a power law fit with apparent fractal dimension Dq ~ 1.93. 
The upper deviation from the power law is at 1/e = 2 10 , i.e. roughly when the number of 
boxes (1/e) 2 is of the order of the number of points N. The dotted and long-dashed lines 
are the power laws with exponents 2 and 1, respectively. 

However, as first remarked in Ref. [11] the slope of the curve in fig. 1 is not constant 
and depends on the resolution. It is already evident in this figure that, as the resolu- 
tion increases (e — » 0), the local slope tends to zero (on this scale most of the points 
are surrounded by very few neighbours). This is better shown in fig. 2a, where the 
discrete approximation of the local slope in logarithmic scale 



is plotted as a function of e for the same set of points. Going from low to high 
resolution, all dimensions from 2 to are covered continuously and smoothly. All 
are apparent fractal dimensions of a set which has indeed dimension 1. It might 
look strange that one does not observe any signal of dimension 1 at any resolution. 
This is due to the inhomogeneity, as shown in fig. 3, where one represents in black 
the occupied boxes at increasing resolutions. While at low resolutions all boxes are 



25 




20 



12 




(9) 



4 



occupied indicating a set of dimension 2, when the resolution increases, the set con- 
centrates on a smaller set of boxes, but already at resolution 1/e = 2 9 one observes 
structures of dimension 1 (lines). However, their weight in the /x-space is negligible 
and, as the resolution is increased even more, the fraction of the set made by spread 
points dominates and consequently dimension 0. A complete multifractal [3] analysis 
should be realized to highlight such inhomogeneities. 

For comparison, we have performed a similar experiment for the Chirikov standard 
map [7] 



pI + i=pI + qI (io) 

4+1 = 4 - Ksm(pl +1 ) mod 2vr , (11) 

where the index t — 0, . . . denotes time and the index i — 1, . . . , NP the number of 
initial particles. Particles are put at t — 0, as above, at q = and with momentum 
uniformly spread in the interval [— p,p}. The phase space is closed by periodizing in 
27r both the q and the p directions and normalized in the unit square. The results 
of numerical experiments conducted above the chaos threshold, at K = 1.2, i.e. a 
parameter region with a mixed (chaotic and regular) phase space, are shown in fig. 2b 
(local slopes) and fig. 4 (occupied boxes in phase space at different resolutions). The 
behavior of all measured quantities is qualitatively the same as the one of the HMF 
model. 




Fig. 2. Local slopes s(e) for resolutions ranging from 1/e = 2 7 up to 2 12 . In Figure a) we 
plot the results for the HMF model at t = 75 with WB initial condition at U = 0.69 and 
N = 10 6 . Figure b) refers to the standard map at t = 160 for K = 1.2 and NP = 10 6 . 

Can one get from straight box-counting an indication of the presence of sets of 
dimension 1? To obtain this results it is enough to perform experiments at lower 
energies for the HMF model and at smaller values of parameter K in the standard 
map. By doing this, one decreases the fraction of phase space which is occupied 
by chaotic orbits and, hence, reduces particle dispersion. In fig. 5 we plot the local 
slopes for the HMF model at U = 0.1 (fig. 5a) and at K = 0.2 for the standard 
map (fig. 5b). In these cases, the local slope, and hence fractal dimension, shows an 
evident plateau at the value 1 as the resolution increases, before dropping to zero 
when the individual points of the set are detected. Let us remark that the number 
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Fig. 3. HMF model. A set of N = 10 points is, initially, uniformly distributed along the 
p axis between — p and p, such that U = p 2 /6 = 0.69. From left to right and from top 
to bottom we show a portion S = [0.33, 0.66[x [0.36, 0.50[ of the //-space, where a box of 
given linear dimension e is colored black if at least one point of the initial condition is 
contained inside it at time t = 75. Resolution ranges from 1/e = 2 7 up to 1/e = 2 12 . 

of particles is the same as before N = NP = 10 6 . This is a strong evidence that the 
line remains a line during time evolution, as expected. 

Having understood that the limited number of particles is crucial in determining 
the results, we show in fig. 6 and fig. 7 the dependence of the local slopes on the 
number of particles for both the HMF model and the standard map. In all cases, s is 
an increasing function of the number of particles, and, at fixed number of particles, 
as the resolution increases, the slope decreases. However, for a more chaotic phase 
space (fig. 6), i.e. HMF at U — 0.69 and standard map at K = 1.2, the curves 
corresponding to the higher resolution show only a slight increase from zero slope 
and one cannot definitely conclude that, as the number of particle increases, they 
will tend to slope 1 (while those which converge, clearly point to dimension 2). On 
the contrary, for a less chaotic phase space (HMF at U = 0.1 and standard map 
at K = 0.2), the curves corresponding to the higher resolutions converge to slope 
(dimension) 1 as the number of particles increases, approximately to (1/e) 2 . 4 



For the standard map we have performed experiments for up to 10 8 particles. We have 
obtained a more evident levelling off of the slope as the number of particles increases and a 
convergence to 1 from above as the resolution subsequently increases. For the HMF model 
the convergence to 1 is followed by a further drop to smaller values at higher resolutions, 
as also shown in fig. 5a: this is again a finite N effect. 
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Fig. 4. Standard map at K = 1.2. A set of NP = 10 6 particles is, initially, uniformly 
distributed along the p axis. From left to right and from top to bottom we show a portion 
S = [0, 0.33[x [0.86, 1.00[ of the phase space, where a box of given linear dimension e is 
colored black if at least one particle of the initial condition is contained inside it at time 
t = 160. Resolution ranges from 1/e = 2 7 up to 1/e = 2 12 . 




Fig. 5. This figure shows how the local slope changes for resolutions ranging from 1/e = 2 8 
up to 2 12 . In Figure a) we plot the results at t = 500 for the HMF model with WB initial 
condition corresponding to U = 0.1 and a number of particles N = 10 6 . Figure b) refers 
to the Standard Map at t = 160 for K = 0.2 and NP = 10 6 . Again the "natural" slope 
s = 1 is roughly at (1/e) 2 ~ N, NP. 

We think that we have solved the puzzle of the existence of a fractal set generated 
in a finite time when starting with a one-dimensional set (a line) in the phase space 
of Hamiltonian systems. The fractal is not there and is a feature of the resolution 
with which the set is observed. 

Of course this does not diminish the role played by fractal concepts in systems with 
long-range interactions. Certainly, when the initial set is already a fractal [12], it's 
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Fig. 6. Local slope as a function of the number of particles (ranging from 10 3 up to 
10 7 ) for different values of the resolution 1/e: 2 7 (A); 2 s (o); 2 9 (x); 2 10 (□); 2 11 (+); 
2 12 (0)- Figure a) shows the results for the HMF model at U = 0.69 and Figure b) the 
corresponding results for the standard map at K = 1.2. 




well possible that the set remains fractal during time evolution. Moreover, when 
studying evolving universe models, one effectively introduces dissipative terms in 
the equations of motion [12,13] and one can indeed observe fractal sets generated 
from non fractals. Another situation is the one of the advection of passive tracers [14] 
or vector fields [15], where one indeed can observe fractal sets in Hamiltonian phase 
spaces. This is due to the presence of non escaping orbits inside the mixing re- 
gion, which form a zero measure fractal repeller [14] or, in the kinematic dynamo 
problem [15], to the concentration of the magnetic flux onto a singular set as t — > oo. 
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